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ABSTRACT 

We compare the structure and kinematics of the 11 known satellites of the Milky Way 
with high resolution simulations of the formation of its dark halo in a ACDM universe. 
In contrast to earlier work, we find excellent agreement. The observed kinematics are 
exactly those predicted for stellar populations with the observed spatial structure 
orbiting within the most massive "satellite" substructures in our simulations. Less 
massive substructures have weaker potential wells than those hosting the observed 
satellites. If there is a halo substructure "problem" , it consists in understanding why 
halo substructures have been so inefficient in making stars. Suggested modifications of 
dark matter properties (for example, self-interacting or warm dark matter) may well 
spoil the good agreement found for standard Cold Dark Matter. 

Key words: galaxies: satellites - galaxies: evolution - galaxies: Milky Way - galaxies: 
dark matter - methods: N-body simulations 



1 INTRODUCTION 

In hierarchical cosmologies such as the now-standard ACDM 
model, objects like the dark halo of our Milky Way grow 
through the merging of previously collapsed systems with a 
wide range of masses. Even the earliest detailed models for 
the growth of a "Milky Way" showed that simple assump- 
tions for star formation efficiency imply many more visible 
satellites than are actually observed (Kauffmann, White & 
Guiderdoni 1993). These authors argued that gas conden- 
sation and thus star formation must be strongly suppressed 
in systems with low escape velocity, perhaps by photoioni- 
sation heating. This echoed Efstathiou's (1992) arguments 
on the related question of why the faint-end slope of the 
observed galaxy luminosity function is much shallower than 
predicted by the simple hierarchical galaxy formation the- 
ory of White & Rees (1978) and its later reworkings in the 
CDM context (e.g. White & Frenk 1991). 

This issue drew relatively little attention until N-body 
techniques became capable of simulating halos with hun- 
dreds of thousands of particles. Moore et al. (1999) and 
Klypin et al. (1999) showed that galaxy halos in a CDM 
cosmogony are not smooth systems. Roughly 10% of the 
mass within their virialised regions is contributed by a host 
of dense self-bound substructures. These are the surviving 
cores of objects which merged together to make the final sys- 
tem, and so correspond directly to the overabundant satel- 
lites of the earlier models. The apparent overabundance may 
again be reconciled with the small number of visible satel- 
lites by invoking the inhibiting effects of photo-heating (Bul- 



lock, Kravtsov & Weinberg 2000, Benson et al. 2002). Both 
Moore et al. (1999) and Klypin et al. (1999) emphasised a 
different problem, however. The maximum circular velocities 
of the 10 or 20 most massive substructures within a Milky 
Way halo are predicted to be in the range 20 to 60 km/s, 
whereas the observed velocity dispersions of 7 of the Milky 
Way's 11 satellites are below 10 km/s. Halo models typi- 
cally have well over a hundred substructures with maximum 
circular velocity above 10 km/s. 

This discrepancy has been considered a "crisis" for 
the conventional ACDM cosmogony, prompting proposals to 
modify the microscopic properties of the dark matter (e.g. 
Spergel & Steinhardt 2000; Moore et al. 2000; Yoshida et al. 
2000; Bode, Ostriker & Turok 2001). We here examine the 
discrepancy more critically, using ACDM simulations with 
similar resolution to the best simulations of Moore et al. 
(1999) or Power et al. (2002). We analyse the potential well 
structure of the most massive "satellites" in our final system, 
calculating the velocity dispersion profiles predicted if stel- 
lar systems identical in structure to the observed satellites 
are placed in each. Remarkably, we find excellent agreement 
with the observed velocity dispersions, even for systems like 
Fornax and Draco where resolved profiles are available (Ma- 
teo 1997; Kleyna et al. 2002). 

As we were completing this work, Hayashi et al. (2002) 
posted a preprint addressing similar issues (and with sim- 
ilar conclusions). Their approach is complementary to the 
one we adopt here. They do not analyse a fully consistent 
simulation of a Milky Way halo and its substructure, nor do 
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they compare in detail with the kinematic structure of the 
observed dwarfs. On the other hand, they carry out a much 
more complete and reliable analysis of the structural evo- 
lution of individual satellite substructures than is possible 
with our own simulation data. We will use their results be- 
low to provide an independent check on the most uncertain 
aspect of our own simulations - whether they give reliable 
estimates for the inner "core" structure of satellite subhalos. 

In the next section we will give a brief description of 
the simulations we have carried out, the methods we have 
used to find substructures and to characterise their potential 
wells, and the checks we have made on the reliability of these 
measures. Section 3 then shows how to predict the line-of- 
sight velocity dispersion profile for a dwarf galaxy of given 
spatial structure within a given potential well. We apply 
the technique to the 11 known satellites of the Milky Way 
and the 20 deepest satellite potential wells of our highest 
resolution simulation. Section 4 discusses our results, their 
principal remaining uncertainties and their implications. 



2 THE iV-BODY SIMULATIONS 

We work with a flat A-dominated Cold Dark Matter uni- 
verse, with matter density Q m = 0.3, cosmological con- 
stant r^A = 0.7, expansion rate Ho = 70 km s _1 Mpc _1 , 
index of the initial fluctuation power spectrum n = 1, and 
present-day fluctuation amplitude as = 0.9. We begin with a 
dark matter simulation of a "typical" region of the Universe 
(N ~ 6 x 10 7 , particle mass ~2x 10 8 M o ) for which the 
techniques of Springel et al. (2001b, hereafter SWTK) have 
been used to follow the formation of the galaxy population 
(Stoehr et al. 2002, in preparation). We identify a relatively 
isolated "Milky Way like" galaxy and resimulate its halo at 
a series of higher resolutions, again using techniques from 
SWTK and the iV-body code GADGET (Springel, Yoshida & 
White 2001a). In the simulations GA0, GA1 and GA2 anal- 
ysed in this paper the resimulated halo has 13603, 123775 
and 1055083 particles respectively within r2oo, the radius 
enclosing a mean density 200 times the critical value. The 
density profiles of these three resimulations agree well and 
are accurately fit by the NFW formula with concentration 
c = 10 (Navarro, Frenk & White 1997). We show the corre- 
sponding circular velocity curves in Figure la. These have 
been scaled down by a factor of 0.91 in velocity and radius 
(corresponding to a factor of 0.74 in mass, but unchanged 
density and time scales) so that they peak at 220 km/s. 
With this scaling, dark matter particle masses are 1.8 x 10 s , 
1.9 x 10 7 and 2.0 x 10 6 M Q , and Plummer equivalent soften- 
ing lengths are 1.8, 1.0 and 0.49 kpc in GA0, GA1 and GA2, 
respectively. In all three V2oo ~ 270 kpc. Note that since the 
Milky Way's stars contribute significantly to its measured 
rotation velocity, our chosen scaling probably produces too 
large a mass for the Milky Way's halo and thus also for sub- 
structures within it. This is conservative for the purposes of 
this paper. 

We identify self-bound substructures within our final 
halos using the procedure SUBFIND described in detail by 
SWTK. In Figure lb we compare the number of substruc- 
tures found within r2oo in each of our simulations. For each 
we plot cumulative abundances down to a mass correspond- 
ing to 20 particles. Clearly the agreement is good. For ex- 



ample, down to 4 x 10 Mq, the 20 particle limit for GA0, we 
find 5, 7 and 2 subhalos in GA0, GA1 and GA2, respectively, 
well within the fluctuations expected for Poisson statistics. 
Down to the 20 particle limit for GA1 (4 x 1O 8 M ) there are 
28 and 27 subalos in GA1 and GA2. We have also verified 
that our cumulative function of peak velocities is consistent 
with that of Font et al. (2001). 

For each subhalo we define the centre as the position 
of the densest particle, where densities are evaluated using 
an SPH smoothing kernel. We then calculate a circular ve- 
locity curve as V c (r) = (GM(r)/r)'' 2 where M(r) is the 
total mass within r of this centre. In almost all of the more 
massive subhalos, these curves rise from the centre, peak 
at a distance of a few kpc, and then fall for a while before 
eventually rising again because of the contribution from the 
smooth halo. Fitting such curves we can define a maximum 
circular velocity and a radius at which it is achieved for each 
subhalo. In Figure lc we plot these two quantities against 
each other for all subhalos in GA1 and GA2 with more than 
100 and 300 particles, respectively. There is a weak correla- 
tion but no clear tendency for the GA1 subhalos to be less 
concentrated (i.e. to have larger r max at given Vm ax ) than 
those in GA2. By comparing the evolution of the two simu- 
lations in detail it is possible to identify counterparts in GA2 
for all the eight GA1 subhalos plotted in Figure lc. Four of 
these lie within r2oo in GA2, and the corresponding pairs of 
points are indicated with ellipses in the plot. There is good 
agreement between the (V m ax,r max ) values measured in the 
two simulations despite the factor of 10 difference in mass 
resolution. 

In Figure Id we plot circular velocity curves for some of 
the 20 GA2 subhalos with the largest values of Umax- Even 
the least massive of these has more than 300 particles, and 
we have checked that they all have maintained their struc- 
ture with relatively little change since at least z = 0.4. It 
is noticeable that these circular velocity curves have nar- 
rower peaks than the NFW form, so we have fitted them 
with parabolae 

log(V c /Vm ax ) = -a [log(r/r max )] 2 , (1) 

adjusting the constant a to get a good fit at r < r max . We 
find values for a ranging from 0.25 to 0.7 with a median of 
0.45. These fits are shown as dotted curves in Figure Id, 
and we use them to represent the circular velocity data in 
the analysis of the following section. This analysis depends 
quite sensitively on the shape of the curves at small radius. 
This is a matter for concern since in this regime, our re- 
sults are likely to be affected by discreteness effects and by 
gravitational softening. Although it is reassuring that we see 
little temporal evolution in the inner structure of individual 
subhalos, and that the more massive GA1 halos agree in 
structure with their GA2 counterparts, a more convincing 
demonstration that equation (1) is realistic comes from the 
fact that it fits the circular velocity curves of the much higher 
resolution simulations of stripped satellites by Hayashi et al. 
(2002). The curves in their Figure 10 are well fit by values of 
a in the range 0.37 to 0.65. Below we will repeat our anal- 
ysis using a = 0.37 for all our substructures, rather than 
the individual values estimated from their detailed internal 
structure. This choice maximises the circular velocity at the 
relevant radii (r ~ 0.1r max ) for curves compatible in shape 
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Figure 1. (a) Circular velocity curves for the "Milky Way" halos in our three resimulations, scaled to peak at 220 km/s. An arrow 
indicates r2oo- We overplot a NFW fit with c = 10 (dotted line). Vertical lines show the Plummcr softening lengths, (b) The cumulative 
number of subhalos within T200 as a function of their mass for these same three resimulations. Vertical lines show the 20 particle limits, 
(c) Maximum of the circular velocity curve against the radius at which it is reached for all subhalos with more than 100 particles in 
GA1 (filled circles) and with more than 300 particles in GA2 (open circles). Vertical lines show the softening lengths and ellipses enclose 
counterpart subhalos. (d) Circular velocity curves and fits to equation (1) for subhalos 1, 2, 5, 8, 11, 14, 17 and 20 in GA2 in order of 
decreasing V max . The vertical line shows the softening length. 



with those of Hayashi et al. and for the values of V ma x and 
r max measured in our simulation. 



3 SUBHALO POTENTIALS AND THE 
OBSERVED MILKY WAY SATELLITES 

Except for the Magellanic Clouds, all the known satellites 
of the Milky Way are dwarf spheroidal galaxies. All have 
measured core and tidal radii, as well as measured central 
line-of-sight velocity dispersions. We summarise the data in 
Table 1. They are taken primarily from the review of Ma- 
teo (1998), but we have preferred the more recent density 
profile parameters measured for Draco by Odenkirchen et 
al. (2001). Sagittarius is also a special case. Here we take 
the core and tidal radii from the pre- disruption model of 
Helmi & White (2001) and carry out our analysis using 
both the current central velocity dispersion given by Ma- 
teo and the higher pre-disruption value from the model of 
Helmi & White. The latter value and quantities based on 
it are indicated in parentheses in Table 1. For Fornax and 
Draco there are also published velocity dispersion profiles 



Table 1. Core radii r c , tidal to core radius ratios rt/r c and cen- 
tral velocity dispersions <ro for the Milky Way dwarf spheroidals. 
The last two columns show the number of simulated subhalos for 
which the predicted central velocity dispersion is larger than that 
observed. The first takes a in equation (1) from a fit to our data, 
the second fixes a = 0.37. 





r c [kpc] 


n/r c 


_ r km 1 

o [— ] 


NGA2 




Sagittarius 


0.44 


6.8 


11.4(19) 


11(2) 


15(1) 


Fornax 


0.46 


5.1 


10.5 


13 


15 


Leo I 


0.215 


3.8 


8.8 


4 


5 


Sculptor 


0.11 


13 


6.6 


4 


5 


Leo II 


0.16 


3 


6.7 


1 


5 


Sextans 


0.335 


9.6 


6.6 


18 


19 


Carina 


0.21 


3.3 


6.8 


6 


10 


Ursa Minor 


0.20 


3.2 


9.3 





2 


Draco 


0.18 


5.2 


9.5 





2 



with good signal-to-noise; we compare these to our ACDM 
predictions below. To characterise the potentials of the Mag- 
ellanic Clouds we take a circular velocity of 50 km/s at a 
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distance of 5 kpc for the LMC (van der Marel et al. 2002) 
and a circular velocity of 60 km/s at a distance of 2.5 kpc 
for the SMC (Stanimirovic 2000). Comparing with Figure 
Id, we see that the Clouds fit well in the two most massive 
subhalos of GA2. We now ask whether the next ten or so 
most massive simulated subhalos could host the observed 
dwarf spheroidals. 

A spherical stellar system with number density profile 
p(r) in equilibrium within a spherical potential well with 
circular velocity curve V c (r) will, if its velocity dispersion 
tensor is everywhere isotropic, satisfy 



dr 



P V C 2 



(2) 



where <r(r) is the one-dimensional velocity dispersion profile. 
Assuming that p(r) goes to zero at the finite "tidal" radius 
r t , this equation can be integrated to obtain 



p(r 



ia 2 (r) 



dr' pVc/r. 



(3) 



Projecting pa and p and taking their ratio, we can obtain 
a p (r p ), the line-of sight velocity dispersion at projected dis- 
tance r p from the satellite centre, 



£<lrpV?(r 



2 /r 



J7 dr pr/(r 2 -r|) 1/2 



(4) 



We apply this formula to stellar systems with the observed 
structure of the dwarf spheroidals embedded in the potential 
wells of our simulated subhalos. We will represent the lat- 
ter by equation (1) and the former by King's (1962) fitting 
formula 



p(r) oc [arccos(z)/z — (1 — z 2 ) 1 ^ 2 ] /z 
where 

l+r 2 /r c 2 



2 

Z = 



l + r t 2 /r 2 



(5) 



(6) 



We have used this method to calculate the expected 
central line-of-sight velocity dispersion of a stellar system 
with the core and tidal radii of each of the observed dwarf 
spheroidals embedded in each of the 20 most massive sub- 
structures in our GA2 halo. For each dwarf spheroidal, we 
count the number of subhalos for which the predicted cen- 
tral dispersion is larger than that observed, and we record 
this number in Table 1. We have carried out this exercise 
twice; once using the fits to the individual simulated cir- 
cular velocity curves (see Figure Id) and once taking r max 
and Knax from the simulation and then using a — 0.37 in 
equation (1) for all subhalos. The latter procedure forces the 
inner structure of each subhalo to resemble that of the most 
concentrated of the stripped satellites simulated by Hayashi 
ct al. (2002). These numbers show a remarkable result. For 
both assumptions all 11 of the Milky Way's satellites can be 
accommodated within the 20 most massive subhalos, and 
all but 3 can be accommodated within the 10 most mas- 
sive (all but 2 if the higher dispersion value is adopted for 
Sagittarius) . Given that substantial realisation to realisation 
fluctuations are expected in the properties of the more mas- 
sive subhalos, there is surprisingly good agreement between 
the kinematics of the observed satellites and those predicted 
by our ACDM simulation. 
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Figure 2. Observed velocity dispersion profiles for Fornax and for 
Draco are compared with those predicted for stellar systems with 
the observed density structure embedded in some of the most 
massive dark matter satellites in our ACDM simulation. Solid 
lines show results for circular velocity curves estimated directly 
from our simulation (GA2), while dotted lines assume a circular 
velocity curve shape similar to that of the most concentrated 
satellites in the simulations of Hayashi ct al. (2002). 



For Fornax and Draco, more detailed kinematic data 
have been published and allow a closer comparison. In Fig- 
ure 2 we reproduce velocity dispersion profiles from Mateo 
(1997) for Fornax, and from Kleyna et al. (2002) for Draco, 
and we overplot the predictions of the above equations for a 
number of our more massive subhalos. From this figure it is 
clear that GA2 not only predicts the correct central velocity 
dispersion values, but also the correct shapes for the disper- 
sion profiles. Within the observational uncertainties this is 
equally true for circular velocity curves taken directly from 
GA2 and for curves with a the most concentrated shape con- 
sistent with the simulations of Hayashi et al. (2002). The 
maximum circular velocities of the subhalos compared with 
the Fornax data range between 31 and 46 km/s. The cor- 
responding range for the subhalos compared with Draco is 
35 to 54 km/s. In a ACDM model the dark halos of dwarf 
spheroidals are predicted to have maximum circular veloc- 
ities much larger than the observed velocity dispersions. It 
is a curious that the observational data require Draco and 
Fornax to have similar halos despite the fact that Fornax is 
more luminous by a factor of almost 60. 



4 DISCUSSION 

The last section shows that the potentials of the most mas- 
sive subhalos in our ACDM simulation are in excellent agree- 
ment with the observed kinematics of the Milky Way's satel- 
lites. There is some reason to be cautious, because subhalo 
concentrations could be significantly underestimated as a 
result of numerical limitations, in particular the relatively 
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large gravitational softening and particle mass in even our 
highest resolution simulation. We have argued from a com- 
parison of simulations with differing resolution that the in- 
duced bias in the values of V ma x and r max is small for objects 
as massive as those which concern us here. Furthermore, our 
subhalo circular velocity curves agree well at small radii with 
those which Hayashi et al. (2002) obtain from much higher 
resolution simulations of the tidal stripping of individual 
satellites. These authors explain why tidal effects produce 
objects with lower central concentration than isolated ha- 
los of similar maximum circular velocity. Even adopting the 
most concentrated profile consistent with their simulations 
has only a minor effect on our analysis. (Compare the rank- 
ings for the two cases in Table 1.) On the basis of the cur- 
rently available observational and simulation data, it seems 
more appropriate to consider the observed kinematics of the 
Milky Way's satellites as a triumph for the ACDM model 
than as a crisis. 

If the observed dwarf spheroidals do indeed have dark 
halos of the kind that both we and Hayashi et al. (2002) sug- 
gest, there are a number of interesting consequences. In the 
first place, the measured "tidal radii" can have nothing to 
do with tidal effects, but must reflect an edge to the visible 
stellar population within a much more extended dark halo. 
This kind of structure appears to be a direct and inevitable 
consequence of the flat or rising velocity dispersion profiles 
measured in Fornax and Draco. Tidal tails and extra-tidal 
stars should not, then, be present in most systems. Although 
there is no problem accommodating a single disrupting ob- 
ject like Sagittarius, it would become uncomfortable if tidal 
stripping were detected unambiguously in other systems. 
There are some indications that Carina might be such a 
case (Majewski et al. 2000). 

A second consequence is that the total mass associated 
with the dwarf spheroidals is much larger than usually as- 
sumed. The tenth most massive subhalo in our GA2 simu- 
lation has a bound mass of 1.2 x 10 9 Mq and its progenitor 
system was several times more massive just before it fell 
into the Milky Way's halo. This sharpens the problem in 
understanding why the dwarf spheroidals have formed stars 
with such low efficiency; their stellar masses lie in the range 
10 6 to 1O 8 M0. Clues presumably lie in their comparatively 
narrow ranges of size and velocity dispersion, and in their 
surprisingly varied star formation histories (Mateo 1998). 

Finally, if the observed satellites do occupy the most 
massive subhalos in a ACDM model, then many smaller sub- 
halos are presumably also present but are devoid of stars. 
It is interesting to look for observable consequences of their 
existence. Dynamical effects, for example heating or distor- 
tion of the Galactic disk or scattering of the orbits of visible 
halo objects, are weak, and are dominated by the most mas- 
sive subhalos rather than by their more abundant low mass 
brethren. Such effects may, nevertheless, be detectable in fa- 
vorable cases (Font et al. 2001; Johnston, Spergel & Haydn 
2001). It may also be possible to detect substructure in grav- 
itational lens galaxies through the statistics of flux ratios in 
samples of multiply imaged quasars (Mao & Schneider 1998; 
Chiba 2002; Dalai & Kochanek 2002); again the effects are 
dominated by the few most massive subhalos. If dark matter 
detection experiments are successful, it may become feasi- 
ble to search for structure in the dark matter distribution 
at the Earth's position. Detailed analysis suggests that in a 



ACDM model the signal will be difficult, but perhaps not 
impossible, to see (Helmi, White & Springel 2002). 

The fact that the core structure of satellite subhalos 
in a ACDM model agrees with the kinematics of observed 
satellites is prima facie evidence against dark matter with 
modified properties, for example, Warm Dark Matter or Self- 
Interacting Dark Matter. In recent work, the primary moti- 
vation for considering such modifications to the microscopic 
physics of the dark matter particles has been to reduce the 
concentration of halos and the abundance of substructure 
(e.g. Spergel & Steinhardt 2000; Yoshida et al. 2000; Bode 
et al. 2001). A significant reduction of the central concentra- 
tion of satellite subhalos in our simulations would, however, 
make it difficult to produce satellites with velocity disper- 
sions as large as those observed. 



REFERENCES 

Bode P., Ostriker J. P., Turok N., 2001, ApJ, 556, 93 
Benson A. J., Frenk C. S., Lacey C. G., Baugh C. M. 

2002 MNRAS, 333, 177 
Bullock J. S., Kravtsov A., Weinberg D 
Chiba M., 2002 ApJ, 565, 17 
Dalai N., Kochanek C, 2002, ApJ, 572, 25 
Efstathiou G., 1992, MNRAS, 256, 43 
Font A. S., Navarro J. F., Stadel J., Quinn T 



Cole S. 



2000, ApJ, 539, 517 



Ha yashi E.. Navarro , : 
astro- ph/0203004 



2001, ApJ, 563, 1 
F., Taylor J. E., Stadel J., Quinn T., 2002, 



Helmi A., White S 
Helmi A., White S 
Johnston K. V 
Kauffmann G.. 



D. M., 2001, MNRAS, 32 3. 529 

D. M., Springel V., 2002, |astro-ph/020128fl| 
, Spergel D. N., Haydn C, 2001, |astro-ph/0111196 
White S., Guiderdoni B., 1993, MNRAS, 264, 201 



W., Gilmore G., Frayn C, 
PradaF., 1999, ApJ, 522, 
E., 



King I., 1962, AJ, 67, 471 
Kleyna J., Wilkinson M. I., Evans N. 

2002, MNRAS, 330, 792 
Klypin A., Kravtsov A., Valenzuela O. 

82 

Majewski S. R., Ostheimer J. C, Patterson R. J., Kunkel W. 

Johnston K. V., Geisler D., 2000, AJ, 120, 2550 
Mao S., Schneider P., 1998, MNRAS, 295, 587 
Mateo M., 1997, in Ariboldi et al. , eds, ASP Conf. Ser., 116, 259 
Mateo M. L., 1998, ARA&A, 36, 435 

Moore B., Ghigna S., Governato F., Lake G., Quinn T., Stadel 

J., Tozzi P., 1999 ApJ, 524, 19 
Moore B., Gelato S., Jenkins A., Pearce F. R., Quilis V., 2000, 

ApJ, 535, 21 

Navarro J. F., Frenk C. S., White S. D. M., 1997, ApJ, 490, 493 
Odcnkirchen M., Grebel E., Harbeck D., Dehnen W., Rix H., 

Newberg H., Yanny B., Holtzman J. et al. 2001, AJ, 122, 

2538 

Power C, Navarro J. F., Jenkins A., Frenk C. S.. White S. D. M., 
Springel V., Stadel J., Quinn T., 2002, astro-ph/0201544 

Spergel D. N., Steinhardt P. J., 2000, Physical Review Letters, 
84, 3760 

Springel V., Yoshida N., White S. D. M., 2001a, New Ast. 6, 79 
Springel V., White S. D. M., Tormen C, Kauffmann G., 2001b, 

MNRAS, 328, 726 
Sta riimirovic S^ 2000. PhD thesis, 

http: / /www. naic.edu/ ~sstanimi/thesis.html , 92 
van dcr Marcl R. P.. Alves . D. R., Hardy, E., Suntzcff, N. B., 



341 



2002, |astro-ph/0205161 
White S. D. M., Rees M. J., 1978, MNRAS, 183, 
White S. D. M., Frenk C. S., 1991, ApJ, 379, 52 
Yoshida N., Springel V., White S., Tormen C, 2000, ApJ, 544, 

87 



© 0000 RAS, MNRAS 000, 000-000 



